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With the assistance of a complex singlet, and an effective operator involving CP violations, the 
dark matter relic abundance and baryon asymmetry of the universe have been addressed simul¬ 
taneously. We studied the electroweak baryogenesis mechanism systematically. The electroweak 
phase transition analysis indicates that the strong first order phase transition takes place by one- 
step or two-step type due to the dynamics of the energy gap between the electroweak vacuum and 
the vacuum of the complex singlet. The relation between the magnitude of baryon asymmetry of 
the universe and the phase transition type and strength has been explored in the framework of 
electroweak baryogenesis. 


I. INTRODUCTION 

The Standard Model (SM) has passed most experi¬ 
mental tests during the last 40 years. And the discovery 
of a 125 GeV SM-like Higgs particle by the ATLAS [1] 
and CMS [5] collaborations at the LHC seems to provide 
the last missing piece of the SM. However, it has long 
been known that the SM has two obvious shortcomings, 
i.e., the explanations of the observed dark matter abun¬ 
dance and matter-anti-matter asymmetry of the universe. 

Firstly, the observed dark matter abundance is one 
issue that couldn’t be addressed in the SM. The exis¬ 
tence of dark matter (DM) has already been established 
by the observations of galaxy rotation curves and analy¬ 
sis of cosmic microwave background (CMB) etc. And the 
PLANCK [3] and WMAP [3] predict that DM constitutes 
about 26.5% of our Universe. Secondly, the matter-anti- 
matter asymmetry of the universe, i.e., baryon asymme¬ 
try in the Universe (BAU) as a baryon to entropy ra¬ 
tio HIS] 

— « (0.7 - 0.9) X 10"^° , (1) 

s 

couldn’t be addressed and predicted in the framework 
of the SM. The dynamical generation of BAU requires 
three necessary ingredients [B]: (1) violation of baryon 
number; (2) violation of both C- and CP-symmetry; and 
(3) departure from equilibrium dynamics or CPT vio¬ 
lation. Though the SM contains all these requirements, 
and the electroweak phase transition (EWPT) provides a 
natural mechanism for baryogenesis [B] , the SM is unable 
to solve the problem for the reasons being listed below: 
Although the baryon number violation is provided in the 
SM by weak sphalerons CHS], the departure from ther¬ 
mal equilibrium, provided by a strong first order phase 
transition, proceeding by bubble nucleation which makes 
electroweak symmetry breaking (EWSB), does not oc¬ 
cur in the SM unj. In the SM, no first order phase 


transition occurs for Higgs mass larger than about 80 
GeV [mil3) . which is far below the experimental bound 
of ruh > 114 GeV from LEP [H] and to^R:: 125 GeV 
from the LHC (TJ [2] . In addition, the CP violation in the 
CKM matrix is too small to produce a sufficiently large 
baryon number. New CP violations beyond the SM are 
required m- 

This work is aimed at accommodating the observed 
DM relic density and the BAU simultaneously. The ob¬ 
served DM relic density reported by CMB anisotropy 
probes suggests weakly interacting massive particle 
(WIMP) j^nmil and a feasible candidate for DM requires 
the extension of the SM. Based on CxSM [52], we con¬ 
sider one simple extension of the SM with a complex sin¬ 
glet (S) being supplemented. The singlet transforms triv¬ 
ially under the SM gauge groups and the imaginary part 
takes the responsibility of being DM. One very attrac¬ 
tive mechanism, i.e., electroweak baryogenesis (EWBG), 
wherein baryon number generation is driven by the CP 
asymmetry at the time of the EWPT^, is adopt. We 
take advantage of the complex singlet S to obtain the 
CP-violation required by EWBG through introducing a 
dimension-6 operator. Ref. |3H] considered the analogous 
dimension-5 operator involving S'/A, and Ref. [32] use 
S^/A^ because of the Z 2 symmetry S ^ —S needed to 
prevent decay of S, as befits a dark matter candidate. In 
our case, we consider the complex scalar S, thus we have 
SS^A^, which ensures that the dark matter candidate 
preserves Z 2 symmetry and the CP violation source sur¬ 
vives after symmetry breaking. The relevant Lagrangian 
takes the form of 

ytQLH (^1 + tR + h.c. (2) 


^ Which should be the strong first-order EWPT to protect the 
generated baryon asymmetry from washout and suppress the 
sphaleron after the process is ended. 
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where a, b are real parameters and A is a new physics 
scale. During the EWPT, the top quark mass gets a 
spatially-varying complex phase along the bubble wall 
profile, which provides the source of CP violation needed 
to generate the baryon asymmetry. Precision tests, es¬ 
pecially the electric dipole moment (EDM) searches, can 
probe directly CP violation relevant to EWBG. Using the 
polar molecule thorium monoxide (ThO), the ACME col¬ 
laboration reported an upper limit on the electron EDM 
(eEDM) recently [16], at 90% confidence level, an order 
of magnitude stronger than the previous best limit, 

|de| < 8.7 X 10-2®ecm . (3) 

This limit severely constrains the allowed magnitude of 
CP-phases in the Higgs couplings [27H30] via Barr-Zee 
diagrams. We would like to mention that, one must take 
into account the tension between the CP-phase being re¬ 
quired to successfully implement EWBG and constraints 
from the eEDM. The preliminary results of the work are 
listed as follows: 

• The EWPT have been explored in two parameter 
spaces: triple couplings (c 2 , ^i) and quartic cou¬ 
plings {d 2 , S 2 ), as well as cubic and quartic cou¬ 
plings (<52, C 2 ). The behaviors of the strength of 
EWPT {vc/Tc) and the energy gap AU ( the en¬ 
ergy difference between the Electroweak vacuum 
and the vacuum of the complex singlet) are found 
to be opposite: a smaller AU corresponds to a 
bigger vdTc- Two types of phase transition, i.e., 
one- and two-step are studied, and the two-step 
phase transition always gives rise to a bigger vdTc-, 
which is found to the same as the real singlet case 
as explored in [23] • The dynamics of the phase 
transition could be characterized by AU to some 
extent. 

• The BAU during the EWPT in the model has been 
explored. The behaviors of BAU as functions of 
triple couplings and quartic couplings match well 
with that of EWPT. The magnitude of BAU is 
proportional to the strength of the strong first or¬ 
der electroweak phase transition (SFOEWPT) for 
the one-step EWPT, and could be a little higher 
in the two-step EWPT case with the sign-flip be¬ 
havior appears at some larger Uc/Tc- 


SFOEWPT situation, and the CP violation phase 
needed for the EWBG is allowed by the constraints 
of ACME. 


II. THE MODEL 

With the notation as in [52], the tree-level potential 
of the model is given by, 

U(i7,S) = ^ 



+ ^ C2^g |g|2 ^ ^ 


where H and S are the SU{2) doublet and complex sin¬ 
glet fields. For simplicity, phases are chosen as: (j)b^ = 0 
(tt), (jiSi = 0 (tt), and 0C2 = 0 (’’■)• To get the min¬ 
imization conditions of the potential, it is convenient 
to represent the SU{2) doublet and complex singlet as 
H = (0, h/d2) and S = (5 -k iA)/d2. Thus Eq. Q 
recasts the form of, 

Vdh,S,A) = + + + 

+ ^ (^2 + ^ 1 ) ^ (^2 ~ ^ 1 ) 

X (52 + A 2 ) + ^52^2 + ^ , ( 5 ) 

0 lb 

in which ( 5 i and C2 can be both positive or negative. As 
could be seen from Eq. (§, the breaking of the global 
U{\) induces one Z 2 symmetry for the imaginary part of 
S, i.e. A, thus makes a feasible DM candidate. Mean¬ 
while, the cubic and quartic interactions that are analo¬ 
gous to [23] are both kept for the EWPT study. 

The three minimization conditions of the potential, 

{dVQ/dh)\h=v,s=x,A=Q = 0, 
{dVQldS)\h=v,s=x,A=o = 0, 
(auo/aA)u^„,s^,,^^o = 0, (6) 

with V = 246.2 GeV and x being the vacuum expectation 
values (VEVs) of the h and S, imply that, 


• The imaginary part of the complex singlet serves 
as the DM candidate. The smaller dark matter 
mass ( ttia) region is severely constrained by the 
direct detection experiment LUX. A larger mA can 
give rise to a relatively higher magnitude of the 
relic density. The EWPT in the higher magnitude 
of the dark matter mass region could be indepen¬ 
dent of ruA- A benchmark scenario which accom¬ 
modates DM and EWBG is presented. In the sce¬ 
nario, the magnitude of the BAU could fall into 
the observed range of PLANGK [3] in the two-step 


2 ^22 *^2 

m =- X - X - V , 

2 2 2 

d2bi 2 / ■\/2 c 2 ^2 2 ^22 

62 = &1 - —V lx - —X - —X - —V 

At the minima, the mass matrix of h and S is, 


|A?;2 \52Xv + 

\52xv + d^v5i \d2x‘^ + d^C2X- 


(7) 


( 8 ) 










3 


and the DM candidate A, with mass = bi — 
^C 2 X — /x, does not mix with h or S. The 

non-zero induces mixing between the SM Higgs 

and the real component of the singlet (S). The mass 
square of the “Higgs-like” particle is given by: 
|(m^+wg+\/( w^ - fo r ml > ml, and 

+ ml- -m2)2 for ^ 2 ^ ^ ^ 

For ml > ml, the “Higgs-like” eigenstate hi will be the 
heavier one, and the “Singlet-like” eigenstate /12 will be 
heavier at the time of ml < ml- Demanding mf^ = 125 
GeV, A could be given as a function of (c 2 , i5i) or {d 2 , S 2 ), 
as shown in Fig. . The mass eigenstates hi,h 2 are 


A 



-600-500-400-300-200-100 0 100 

(5i(GeV) 

A 



FIG. 1: Top: A as a function of Si and C 2 from = 125 GeV 
with relevant parameters being set as: x = 200 GeV, rriA = 
300 GeV, d 2 = 1.4, S 2 = 0.2; Bottom: A as a function of S 2 
and d 2 with rrihi = 125GeV, relevant parameters are set as: 
X = 200 GeV, niA = 300 GeV, = —350 GeV, and C 2 = 100 
GeV. 


related with gauge eigenstates h and S as, 


hi 


coscj) 

sin^ 


h 

.^ 2 , 


— sin(() 

cos (j) 


S 


( 9 ) 


The mass eigenstates hi, /i 2 couple to the fermions and 
gauge bosons via SM Higgs couplings reduced by a factor 
of cos (j) and — sin cj) respectively. We take the mixing 
angle to be —7r/4 < 4> < 7r/4 so that hi{h 2 ) is always 
the “Higgs-like” (“singlet-like”) eigenstate. Using this 
convention, the mass matrix elements of Eq. ([^ could be 
expressed as, 


ml = cos 4>'^ml^ + sin (fy'^ml^ , (10) 

TOg =sm(p mf^_^ +cos(p mj^^ , (11) 

mis = cos ^ sin J , (12) 

which allows us to replace A, d 2 and di with mhi,mh 2 
and (f>. As for the vacuum stability condition, parameters 
living in (i 2 < 0 and 62 < —i/\d2 [221 [24] will be excluded. 
For more studies on the stability problem in this kind of 
model, we refer to [2iH25] . 


III. EWBG ANALYSIS 


Before the EWBG is triggered, the early Universe is 
assumed to be hot and radiation-dominated, containing 
zero net baryon charge, and the SU{2)l x U{1)y elec- 
troweak symmetry is manifest |33H36j . With the uni¬ 
verse cooling to temperature below the electroweak scale, 
the Higgs field acquires a non-zero expectation value and 
thus spontaneously breaks the SU[2)lxU{\)y symmetry 
down to U{1)em- During this phase transition process, 
the EWBG takes place. 

Bubble (of the broken phase) formation and growth 
begin at nucleation temperature Tjv( < Tc, the critical 
temperature of SFOEWPT). Baryon creation in EWBG 
takes place in the front (and vicinity) of the expanding 
bubble walls m- With CP violation implemented, par¬ 
ticles in the plasma (symmetric phase) scatter with the 
bubble walls and generate CP (and C) asymmetries in 
particle number densities; and these asymmetries diffuse 
into the symmetric phase ahead of the bubble wall and 
they bias electroweak sphaleron transitions m to pro¬ 
duce more baryons than antibaryons; and some of the 
net baryon charge being created is then captured by the 
expanding wall into the broken phase; at last, inside the 
bubbles(broken phase), the rate of sphaleron transitions, 
which provides (B -|- L)-violating processes and could 
wash out the baryons created, should be strongly sup¬ 
pressed, i.e., the SFOEWPT is needed. 
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A. SFOEWPT realization 


The dynamic of the EWPT is governed by the effec¬ 
tive potential at finite temperature, i.e., 14// (h, S, A, T) 
in our model. After minimizing Veff{h,S,A,T), the 
VEVs of three scalar fields h, S, and A at each tem¬ 
perature could be obtained. When a first order phase 
transition takes place, a local minimum with {h) ^ 0 
develops and becomes degenerate with the electroweak 
symmetric minimum as the temperature decreases to the 
critical temperature Tj,, and the two minima are sepa¬ 
rated by an energy barrier. To protect the BAU gener¬ 
ated by EWBG from being washed out, the strong first 
order condition 


> 1 , (13) 

with Uc being the Higgs VEV of the symmetry broken 
minimum, is necessary. Here, when one uses the condi¬ 
tion, the following three most important theoretical un¬ 
certainties need to be kept in mind [18j : Uc in the condi¬ 
tion is indeed the value of the classic field at minima at 
Tc{(j){Tc)) which depends on the gauge choice; the gauge 
dependence of Tc does not compensate that of (j){Tc)\ and 
the choice of the washout factor is subject to additional 
uncertainties. 

The standard analysis of Veff{h,S,A,T) involves 
tree level scalar potential {Vo{h,S,A)) and Coleman- 
Weinberg one-loop effective potential {Vi{h,S,A)) at 
T = 0, the counter terms Vet being chosen to main¬ 
tain the tree level relations of the parameters in the 
Vb(h, 5', A), and the leading thermal corrections being 
denoted by V (h, S, A, T), 


He// (h, S, A, T) = Ho {h, S, A) + H (h, S,A) + V {h, S, A, T) 
+ Vetih,S,A) . (14) 

Here, Hq is the same as the tree level potential Eq. ([^, 

Hi is the Coleman-Weinberg one-loop effective potential 
at T = 0 and could be expressed in terms of the field- 
dependent masses rrii {h, s, A): 


jr ti. c rnf(h,S,A) 

Vi(h,S,A) = }^n, - ^ - 


In I j - C. 


_ (15) 

in MS scheme and Landau gauge. Here, the sum i 
runs over the scalar, fermion, and vector boson contri¬ 
butions. rii is the number of degrees of freedom for the 
f-th particle, with a minus sign for fermionic particles. 
Q is a renormalization scale which we fix to Q = v. 
We take Ci = ^ {Ci = |) for the transverse (longitu¬ 
dinal) polarizations of gauge bosons, and Q = | for all 
other particles. The field-dependent masses are given in 
the Appendix. B. With Hl being entailed in the poten¬ 
tial, the minimization conditions Eq. will be shifted 
slightly, and the relations Eq. Q do not hold as well. To 


maintain these relations, the counter terms Ht should be 
added [38l [39] , and cast the form of. 

Vet = Smfh^ + 5m\S'^ + 5m\A^ , (16) 


in which, the relevant coefficients are determined by. 


5m\ 

5m\ 

5m\ 


1 dVi 

2v dh , 

(17) 

1 dVi 

2x ds , 

(18) 

1 a^Hi 

2 \h—v,S—x,A—0 • 

(19) 


Thus, the VEVs of h and S as well as the dark matter 
mass will not be shifted. We do not include more com¬ 
plicate terms to compensate the shift of mass matrix of 
S and h, because these shift effects are basically small. 

The finite temperature component V{h,S,A,T) is 
obtained as 


V{h,S,A,T) = ^Y.^MB 


27r^ 


f Mi{h,A,S,T) 
I f2 


( 20 ) 


where the functions JB,F{y^) are given by 


JB,F{y) = j dxx^ \n 1 =F exp (^-y/x"^ + y^ , ( 21 ) 


with the upper (lower) sign corresponds to bosonic 
(fermionic) contributions. And M^{h, S, A,T) is given 
in terms of the field dependent mass at T = 0 and the 
finite temperature mass function H^ as. 


Mlih,S,A,T) = mUh,S,A) + nk , ( 22 ) 


with Hfc being listed in Appendix. B. 

The behaviour of EWPT strongly depends on the 
shape of the potential at zero temperature. Eig.[^ shows 
contours of Vq in the h — S plane in two typical situa¬ 
tions. Since the Vq does not develop a minimum at A 7 ^ 0 
when the dark matter mass is large as will be explored in 
Sec. jlVB] we conduct the analysis in the A = 0 plane. In 
the top panel, the electroweak vacuum is the only local 
minimum of the potential Vq. For this situation, the uni¬ 
verse might transform from the high-T symmetric phase 
to the electroweak symmetry broken phase directly. In 
the bottom panel, an additional local minimum besides 
the electroweak vacuum appears in the direction of S at 
h = 0. For this situation, there is a possibility that the 
universe transforms to this phase at the first step and 
then to the electroweak symmetry broken phase at the 
second step. And in fact, the first investigation of the 
real singlet case [23] shows that the two-step transition 
(with the second step from the high temperature singlet 
vacum to the electroweak vacuum) gives rise to a higher 
magnitude of EWPT strength in comparison with the 
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FIG. 2: The contours of Vo (in units of GeV*) in h — S ( 
in units of GeV ) plane. With parameters being fixed as: 
X = 200 GeV, 5i=-300 GeV, d2=1.4, <52=0.2, A=1.0, and 
C 2 = —10 (100) GeV for the top (bottom) panel. 


one-step transition case. ^ Since the potential at this 
additional minimum (in the direction of S) can be very 
close to that of the electroweak vacuum, the second step 
can happen at a very low temperature, therefore a strong 
first order phase transition could be realized easily. Of 
course, the phase transition could also be the one-step 
case, as long as the energy gap between the minima in 
the direction of S and h is too large to be overcome by 
the T dependent thermal correction at high temperature. 
For the similar analysis in a triplet scenario and Next-to- 
Minimal Supersymmetric Standard Model (NMSSM) we 


^ It is also feasible to study fermion dark matter together with 
strong first order phase transition m- 


refer to [H] and [421144) . 
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FIG. 3: EWPT with parameters being fixed as: x = 
200 GeV, niA = 300 GeV, S 2 = 0.2, d 2 = 1.4, and A be¬ 
ing solved from rrihi = 125 GeV as shown in Fig.[^ The blue 
shade is the regions where ruhi = 125 GeV couldn’t be ob¬ 
tained. The grey region depicts that the electroweak vacuum 
is a local but not global minimum. The light orange and red 
regions are parameters space where one-step and two-step first 
order phase transitions take place. The second order phase 
transition happens in the light blue region. Top: The green 
line is to divide the parameter spaces into the regions where 
the minimum in the direction of S exists (above) and disap¬ 
pears (below). The black contours in the region above the 
green line are the contours of AV = Vs — 14, in units of 10® 
GeV^; Bottom: The contours of Uc/Tc, coscf), and are 
shown with orange(and red), pink and brown lines. The re¬ 
gion excluded by cos (j> < 0.8, which characterizes the mixing 
between h and S, is covered by light grey shade. 

We take advantage of the CosmoTrainsitions pack¬ 
age |45j to investigate the phase transition numerically. 
Fig. [3] shows the result in the C 2 — (5i plane. In the top 
panel, the contour lines of Vq are plotted to help under¬ 
stand this pattern qualitatively: for parameter regions 
bellow the green line, the electroweak vacuum is the only 
minimum of Vq, thus only one-step phase transition could 
happen. The contour AV = 0.6 x 10® GeV'* separates the 
one-step from the two-step type phase transition. When 
the difference between the two minima is bigger than 
that value, the T dependent thermal corrections could 
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only changes the magnitude of ISV a bit thus only the 
one-step phase transition could be obtained. When the 
magnitude of AV^ gets smaller, the phase transition type 
changes from one-step to two-step. The bottom panel 
indicates that, with the increase of C 2 and 5i, one gets 
smaller VcjTc for the one-step phase transition, and a 
second order phase transition is obtained at last. For the 
two-step cases, vdT^ gets larger with the increase of C 2 
and i5i and the SFOEWPT condition Uc/Tc > 1 could 
always be satisfied. The pattern of the contours of VcjTc 
is opposite to that of AE, i.e., the smaller AV leads to 
the larger Vc/Tc- 




FIG. 4: With parameters being fixed as a; = 200 GeV, ttia = 
300 GeV, (5i = —350 GeV, and C 2 = 100 GeV. The descrip¬ 
tions of the two pictures are the same as in Fig. 


LHC |3I], while most two-step ones survive. And to 
make a decisive conclusion on this topic, more data are 
expected from the run II of the LHC as well as other 
future colliders, such as the International Linear Col¬ 
lider (ILC), TLEP(renamed as FCC-ee recently), China 
Electron Positron Collider (CEPC), or the Super proton- 
proton Collider (SPPC), etc. Here, we would like to men¬ 
tion that the electroweak precision observables are sup¬ 
posed to constrain the Higgs and singlet mixing angle (j) 
and the heavy Higgs mass mh 2 [22H241I67], which would 
also bound the SEOEWPT and BAU favored regions. We 
plot the contours of the mixing angle and masses of the 
“Singlet-like” particle in our parameter regions to take 
into account this thing. For detailed analysis we refer 

to [SI]. 


B. BAU gereratation 


To get a more comprehensive understanding of 
EWBG, we solve the transport equations as explored 
in Ref. [ID] directly, in which the speed and profile of 
the walls are the core ingredients. With the nucleating 
bubbles reaching a sizable extent and expanding with a 
constant velocity, it is appropriate to boost into the rest 
frame of the bubble wall and proceed all calculations with 
a planar wall. With z being the coordinate transverse to 
the wall, h{z), S{z) and A{z) are the expectation values 
of h, S and A inside and outside the bubbles. To deter¬ 
mine them precisely, one needs to minimize the Euclidean 
action |iS] : 


S' = 


dz 


' —OO 
1 


+ \{d^Sf 


A^{d.AY + Veff{KS,A,Td 


(23) 


As a good approximation, we estimate the bubble wall 
profiles in the form of. 


KY = - tanh(z/L„)) , (24) 

S{z) = Xc + ]^Axc{l + tajah{z/Ld) , (25) 

A{z) = Oc-I- ^Aac(l -I-tanh(z/Lu,)) . (26) 


Fig. ID illustrates the dependence of the behavior of 
phase transition on the quartic couplings 62 and d 2 - For 
the one-step phase transition, the magnitude of Uc/Tj, 
becomes higher and higher with the decrease (increase) 
of 82 (^ 2 ). For the two-step type phase transition, the 
value of Vc/Tc increases with the two couplings getting 
larger. 

As depicted by the bottom panels of Fig.j^and Fig.|^ 
for our parameter choice, one-step type phase transition 
parameter regions are severely constrained by the mix¬ 
ing angle by fitting current Higgs measurements at the 


Here, is the width of the wall, Vc, Xc, and Oc are 
the VEVs of the electroweak symmetry broken phase at 
Tc, and Ax^ (Aoc) is the total change of (S') ((A)). For 
the scenario of this work, A hardly gets VEV during the 
whole phase transition process, therefore, A(z) is zero, as 
will be explored more in sec tion |IVB[ The wall wid th is 
approximated with « d + ^c) /814 P?] - 
in which Vx is the height of the potential barrier at Tc. 
The bubble wall profiles are shown in Fig. 

As mentioned in Sec. |T| we introduced the CP violat¬ 
ing top quark mass term, and for definiteness we fix a = 0 
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new phase old phase 



FIG. 5: The wall profiles at Tc with new (old) phase corre¬ 
sponding to Electroweak broken (unbroken) phase. 


to maximize the CP violation. The top quark acquires 
a spatially varying complex mass inside the bubble walls 
during the phase transition, being obtained as 

■mtiz) = ■ 

(27) 

following the method used in. [48l |49j. We will verify if 
the nontrivial phase 9{z) is sufficient to source the baryon 
asymmetry in the following studies. The energy scale will 
be chosen as A = 1000 GeV in relevant calculations of 
this work, i.e., eEDM and BAU.^ 

With the top quark profile of our model being ob¬ 
tained, we solve the transport equations ( see Apendix 
D for more details) to get the chemical potentials of 
left-handed SU(2) doublet top quark( fit, 2 ), left-handed 
SU(2) doublet bottoms {fib, 2 ), left-handed SU(2) singlet 
top quark {fit’=, 2 ), and Higgs bosons {fih, 2 ), and the cor¬ 
responding plasma velocities. The chemical potential of 
left-handed quarks cast the form of, 

MSi ~ 2 ^Kx,t)flt,2 + 2 ^Kx^b)flb,2 — 2A'i_(/iic_2 , 

(28) 

where K factors are thermal averages [40]. 

Finally, The baryon asymmetry is obtained using 



405r^, 

‘in'^v^g^T 


dz fiBL{z)e 


(29) 


where P^js is the weak sphaleron rate and v = 
A'bTu]s/{^Vw)- The effective number of degrees of free¬ 
dom in the plasma is g* = 106.75. And we take the wall 
velocity to be =0.1 in our calculation. 


® Here, we would like to mention that the energy scale A charac- 
terizes the size of 0{z) to some extent. And a larger A always 
lead to a smaller 0{z) and results in a lower magnitude of rfs 
and de/e. The dimension six operator in Eq. § induces negligi¬ 
ble effects in the EWPT process, and the wall profiles shown in 
Fig. [^ensures that the operator’s contribution to mt{z) is small 
in comparison with the renormalizable one m- 
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FIG. 6: The correlation between baryon asymmetry and the 
EWPT parameters. Parameters are fixed as Vc = 150 GeV, 
ah = 0, and Tc = 100 GeV(xc = 150 GeV, Xh = 0, and 
ah = 100 GeV) for the top (middle) panel. The magnitude 
of gB as a function of Xc and Xh, Vc and Tc, and VcjTc are 
shown in top, middle, and bottom panels . 





















To get a first glimpse of the behavior of the baryon 
asymmetry with respect to the parameters of EWPT ( 
i.e., the critical temperature and the field expectation 
values Vc,Xc,Xh---), we perform the calculation of BAU 
with these parameters being input directly. We need to 
note that the values of Vc and Tc used here are not phys¬ 
ical parameters we get during EWPT as explored in the 
previous section. The top panel in Fig. [^illustrates how 
the expectation values of S affect the magnitude of BAU, 
with Xc and Xh denoting the VEV of S at electroweak 
symmetry broken phase and symmetric phase. It depicts 
that r]B gets larger with Xh getting smaller for Xc > 0, 
and the situation changes to be opposite when Xc < 0. 
rjB changes its sign near Xh = Xc- In compare with that 
of Xh, r]B relies on Xc more, especially when \xc\ is small. 
The middle panel indicates how r]B depends on Vc and 
Tc. The dependence of riB on Vc is simple and obvious: a 
larger value of riB corresponds to a larger Vc as expected, 
since the source of CP violation depends partly upon the 
amount of variation of |mt|(~ Vc) inside the bubble wall. 
The dependence oi rjB on Tc is more complicated, ijb 
have a almost linear correlation with Tc for larger Tc, 
and the sign oi rjB flips for smaller Tc. The bottom panel 
gives the magnitude of BAU as a function of fc/Tc, which 
indicates that some value of Vc might induce unphysical 
results, i.e., rjB < 0. Thus, a larger v^Tf, may not leads 
to a bigger rjB- 


Then, we use the results obtained in the section, [ill A| 
as the input parameters and calculate the magnitude of 
the BAU being generated during SFOEWPT. The results 
are shown in Fig. which depicts that for the one-step 
case the behaviors oi rjB are very regular and indicates 
one good correspondence with the behavior of Vc/Tc, i.e., 
a simple linear correlation. While, a more complicated 
correlation for the two-step case exists. The sign of tjb 
flips and flips again with the increase of VcjTc in the 
two-step case, which depicts a similar behavior as that of 
Fig.jej 

As could be observed in the two panels of Fig. a 
larger baryon asymmetry can be generated during the 
two-step type phase transition process. And, 7 x 10“^^ < 

< 9 X 10“^^ turns out to live in the two-step type 
EWPT favored region for the parameters set. Last but 
no least, we would like to go more deeper to the inner 
mechanism of the sign flip oi rjB and relevant behaviors 
with respect to Vc/Tc as shown in Fig. [^ The BAU mag¬ 
nitude is highly related with that of jiB^ (see Eq. 291, 
which is supposed to be dominated by the behavior of 
the source term St ■ To illustrate the physical picture 
more clear, we take some typical benchmark points be¬ 
ing used in Fig. [^ and plot the magnitude of and 
the dominated source term St with respect to the spa- 
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FIG. 7: The BAU magnitude rjB as a function of VcjTc, with 
orange(red) points depicting the one-step (two-step) phase 
transition cases. Parameters are set as: x = 200 GeV, rriA = 
300 GeV, 52 = 0.2 ( 5i = -350 GeV) and da = 1.4 
(ca = 100 GeV) for top (bottom) panel. . 


cial coordinate 2 , see Fig.[^ The oscillation amplitude of 
jjLBi^ and St are found to be increasing with the increase 
of VcjTc. With the increase of the oscillation amplitude 
of jiBc ^-Tid St, rjB increases at first and decreases lat¬ 
ter. And the volatility of the oscillation of jiBc ^-nd St 
give rise to unphysical BAU magnitude, i.e., the negative 
value of rjB appears (see the red line in Fig. [^. 


C. EDM 


^ Which is the same as that of the Sq in HQ], here we denote it as 
St to depict that it is the top source term. 


In this section, we use the eEDM search to constrain 
the CP violation phase. The eEDM contribution is dom¬ 
inated by the Barr-Zee diagram [^, as illustrated in 
Fig.§ Where, the two scalar mass eigenstates(/ii^ 2 ) give 
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Z 


<52=0.3,£;2 = l-8,Vc=235.3,rc=81-3,;;B=-15.3xlO“" 
i52=0.4,d2 = 1 .2 ,Vc=234.2,7’c=84.3,?;b=0.04x 10“'' 
52=0.3,<;2 = 1.4,vc=222.7,r£, = 101.7,;;5=10.0xl0“" 

<52=0.4,rf 2 =0.8,vv=215.78,7’c=m.9,i;fi = 18.7xl0“' 
<52=0.2,</2 = 1.2,Vc=179.6,r£, = 136.1,<75=9.77x10“" 


FIG. 8: Top and middle panels are the plots of and St 
as a function of the coordinate transversing to the wall (z), 
with benchmark points being given in the bottom panel. 


where TOe,t are the electron (top quark) masses, and Z 
characterizes the mixing between h and S |52j , being ex¬ 
pressed as: 


Z = cos d> sin 6 = 


a /{ml — m|)2 -I- 4 to| 


(31) 


Sh 


and the function G(z), 
<•1 


G{z) = zf 
Jo 


dx- 


■(l-x)-. 


log 


:{1 — x) 


(32) 


The magnitude of eEDM, in the parameter regions of 
Fig. decreases with the increase of (5i((52) and is below 
the limit given by the experimental results of ACME m, 
which means that the CP violation phase in the param¬ 
eter spaces is allowed by the eEDM experiment ACME. 


IV. DARK MATTER AND ELECTROWEAK 
BARYOGENESIS 


In this section, we investigate the possibility to ex¬ 
plain the dark matter relic abundance and EWBG in the 
same parameter spaces in the model. To achieve this 
purpose, we go through the dark matter issues at first. 
Then we explore the relation between the dark matter 
mass and the EWPT. The explanation of the the dark 
matter relic abundance and the EWBG will be acquired 
simultaneously at last. 


two opposite sign contributions, obtained using ® 


A. Dark matter 


de 

e 


16 a TTle 

3 (47r)3 


X 



G{rriflml) 


G {nitlml) 


(30) 


^ Here, only the dominate top quark contribution and ^77 inter¬ 
action have been considered. As for the most general case, and 
possible cancellation mechanism we refer to ED- 


In our model, the thermal relic density flAh^ is mainly 
controlled by the pair annihilation cross section of A, i.e., 
ilAh^ ~ ^/{o'v), in which (av) is the thermal average 
of the product of annihilation cross section and relative 
velocity, with relevant Feynman diagrams being shown 
in Fig. 

Direct detection experiments are searching for dark 
matter scattering off atomic nuclei, and in the model we 
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FIG. 10: Annihilation channels of dark matter particle pairs. 


only have the Higgs mediated spin-independent scatter¬ 
ing cross section which constrains our parameter spaces. 
The scattering cross section of A with a proton for our 
model casts the form of, 


CTdd = 


gAAhi cos ( 


9AAh2 sin 1 


2Trv'^ {iTip + niA) 




X ( fpu + fpd + fps + (3/g) 


(33) 


which is the same with Eq. (25) in [21] except that the 
couplings gAAhi nnd gAAh 2 are different as shown in the 
Appendix. B. 

Fig.[TT] depicts the relic density as a function of tua 
and 62 (C 2 ). The magnitude of the relic density in¬ 
creases with the increase of since the annihilation 
cross section is inversely proportional to the dark matter 
mass. The LUX experimental results mostly constrain 
the smaller dark matter masses region, because the add 
has the dark matter mass suppress effect, as could be 
seen in Eq. [^ The 62 and C 2 affect the magnitude of 
the relic density and direct detection of the dark matter 
mainly through the couplings gAAhi ^-ncl gAAh 2 ■ 


B. The dark matter mass and the EWPT 





rriA (GeV) 



iha (GeV) 


FIG. 11: Dark matter relic density as a function of C 2 

[ 82 ) and niA with x = 200 GeV, (j) = 0.6, and mh 2 = 250 GeV. 
The contours of Q,Ah? are shown by purple lines with C 2 = 30 
GeV {82 = 0.2) for the top (bottom) panel. The parameter 
spaces excluded by LUX [60] are shown by the orange region. 


In this section, we explore the relation between the 
dark matter mass and the EWPT. 

• The dark matter mass, T = 0 and T 7 ^ 0 vacuum 
structure 

The full analysis of the vacuum structure in our 
model is not easy, since there are many extreme 
points in the three dimensional space and the lack 
of Z 2 symmetry of S makes these extreme solutions 
very complicated. In the numerical scan of the 
parameter space in the EWPT part, we find that 
in most cases there are no local minimum at A 7 ^ 0 , 
whether at T = 0 or T 7 ^ 0 , especially when tua 
is large. This can be understood partly through 


the following arguments. Suppose there is a local 
minimum (vh,Vs,VA), it should also be a minimum 
in the S = Vs,h = Vh subspace. In this subspace, 
the potential Vo(vh,Vs, A) has a local minimum 
when the coefficient of A^ term is negative: 


62 + 61 V 2 c 2 . c( 2 2 I 62 2 , n 

+ jVd<0. (34) 

Using the Eq. [^ we rewire the above condition as 


m^ + 


\/ 2 c 2 


g (a:-Us) + —(v.-x ) + — (vd-v ) < 0 . (35) 


As a nature expectation, Vg ^ x and Vh ~ v. So 
when ruA is large, the condition of Eq. 35 can be 
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easily violated and a local minimum at ^ 7 ^ 0 does 
not exist in this subspace, neither in the whole 
space. At T 7 ^ 0, thermal contributions stabilize 
A field at the origin so A 7 ^ 0 minimum does not 
exist at high temperatures. 

• EWPT and the dark matter mass 


could be obtained in some parameter spaces of S 2 and 
C 2 . As for the strength of EWPT and the magnitude 
of BAU, we found that the two-step SFOEWPT case is 
more favored. The magnitude of eEDM is obtained to 
be |(ie| = 4 X 10“^® ecm, below the upper limit given by 
ACME. 


With A = 0, the Eq. [precasts the form of 

Vo(h,S) ="fh^ + ^h^ + ^Sih^S + ^h^S^ 

4 16 8 8 

+ ^(b2-b^)S^ + ^C2S^ + ^S\ (36) 

Thus, the vacua gap AV = Vh — Vs at tree level 
could be obtained as. 


AP = - 


rrVv'^ 1 


4 4^ ^ ® 6^2 16 


b2 — hi 0 v^x5i 

-. ^ -^ 

4 4v^ 


v‘^x‘^S2 A 


- V 

16 


(37) 



62 


with Vs being the solution of 95Po(0, S) = 0. Since 
the dark matter mass could only affect the com¬ 
bination bi -|- 62 , and the combination 62 ~ bi is 
a constant with other parameters being fixed (see 
Eq. , one finds that the dark matter mass is ir¬ 
relevant to the vacuum energy gap, thus one could 
expect that the strength of EWPT (inversely pro¬ 
portional to AP) won’t be affected by which. 


FIG. 12: The contours of the relic density of the dark matter 
{Q.Ah’^) and the magnitude of Vc/Tc with x = 200 GeV, (f) = 
0.6, niA — 1500 GeV, and = 250 GeV. The contours of 
Q,Ah^ are shown by purple lines and the light orange(red) lines 
depict the value of vdTc- The green points are parameter 
spaces with appropriate BAU magnitude 0.7 x 10“^° < rjB < 
0.9 X 10~^°. 


Based on the above two arguments, we found that the 
EWPT could be independent of the dark matter mass in 
the relatively larger dark matter mass regions. 


V. COMMENTS ON THE CALCULATION 
METHODS INVOLVING SFOEWPT AND BAU 


C. EWBG and DM 

The higher magnitude of the relic density requires 
a relatively larger dark matter mass, about tua 1000 
GeV. And in this mass region, as explored in the last sec¬ 
tion and the section [IlI B[ the EWPT and the generation 
of BAU is independent of ttia- In order to accommodate 
the dark matter relic abundance and the EWBG simul¬ 
taneously we assume the dark matter mass tua = 1500 
GeV. 

We present in Fig. [T^the combined result of EWBG 
and DM relic density. All parameter regions in this 
plot are under the direct detection constraint from LUX. 
Fig. [^indicates that the behaviors of one- and two-step 
phase transitions are dominated by S 2 and C 2 respec¬ 
tively. The dark matter relic density could be affected 
by S 2 and C 2 , since they enter into the couplings of the 
dark matter particle and Higgs, i.e., gAAhi{h 2 )i which de¬ 
termines the dark matter particle pairs annihilation pro¬ 
cesses. A larger value of C 2 leads to a larger VLaH^, and 
the correct dark matter relic density observed by mi 


As pointed out by [18], the SFOEWPT criteria Vc/Tc 
calculated in the traditional mathod is subject to gauge 
problems, as the Vc and Tc both involve gauge depen¬ 
dence problem. The EWPT is performed in Landau 
gauge in this paper and one needs to take this into ac¬ 
count when interpreting our results. We have checked the 
robustness of the analysis in the model by following the 
gauge independent treatment introduced in [18] for part 
of the parameter spaces favored by the EWBG and DM 
explanations in Fig. We find that the magnitude of 
Tc drops a lot and the corresponding u(Tc)/Tc is higher 
than Uc/Tc for the two-step case. 

The numerical calculation of the BAU is based on the 
source terms and chemical potentials transport equations 
being computed within the WKB approximation. The 
results being given might be different to some extent in 
comparison with the results calculated using the quan¬ 
tum transport equations from the close time path (CTP) 
formulation, wherein the GP violation source is derived 
using the “(vacuum expectation value)vev-insertion” ap¬ 
proach [71^. 
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VI. DISCUSSIONS AND CONCLUSION with 


By supplementing a complex singlet to the SM, DM 
relic abundance and EWBG could be accommodated si¬ 
multaneously. Both two ingredients of EWBG, i.e., the 
SFOEWPT and new GP violation phase, have been ex¬ 
plored at the same time. The imaginary part of the com¬ 
plex singlet (A) makes the DM candidate due to the re¬ 
duced Z 2 symmetry from a global 17(1) symmetry. 

The phase transition could be one- and two- step 
types as found by the study of the energy gap AV. In or¬ 
der to reveal the effects of cubic and quartic terms on the 
EWPT, we performed EWPT analysis numerically using 
CosmoTransitions package in the parameter spaces of 
(c 2 , ^ 1 ), (^ 2 ) < 52 ) and (^ 2 , C 2 ). Both one- and two- steps 
first order phase transition are obtained, and the two- 
step case favors SFOEWPT more than that of one-step 
case. And this study is consistent with the analysis of 
AV, i.e., a lower magnitude of AV results in a larger 
strength of EWPT. 

The newly introduced GP violation makes the EWBG 
mechanism viable, and the CP violation phase is found 
to be allowed by constraints of ACME. The general and 
specific calculations oi rjs shows that the magnitude of 
BAU is highly related with Vc and T^, and the increase of 
VcjTc does not always cause the increase of tjb, because 
Vc/Tc affects the CP violation source. 

The magnitude of the relic density of the DM is found 
to increase with the increase of the dark matter mass. 
The smaller DM mass regions are severely constrained by 
direct detection experimental results. And the EWPT is 
irrelevant to the DM mass in the larger DM mass regions. 
At last, we present a scenario in which both correct DM 
relic density and successful EWBG could be achieved 
in some parameter spaces of <52 and C 2 considering the 
eEDM experiment ACME. 

There is a possibility that baryogenesis and dark mat¬ 
ter are generated at the time of EWPT simultaneously 
with the dark matter sector being asymmetric dark mat¬ 
ter pMSH] , thus shed light on the path to explain the ra¬ 
tio between the observed abundance of DM and baryons 
in the present day: Qom/^b « 5.4 [3]. The key in¬ 
gredient for the test of this kind of model is the colliders 
search of the Higgs triple coupling [STJ [MU591IHT] and the 
introduced CP violation [62ll66j . that will be postponed 
to the future works. 


Appendix A: Field dependent masses 


= 


mlhih,S,A) mls{h,S,A) ml^{h,S,A) 


Tnlsih^S^A) mlg{h,S,A) ml^{h,S,A)) 

i^hA -S': rrisA -S': "m-AA -S', A) 


and matrix entries are given by 


'-hh 


'^\s 






^‘sA 


rn^A 


m 

v^<5i 


7|lls + |(s“ + yl=) + ?lftqA 2 ) 


4 


-hi + 62 




A- 


•\/ 2 C 2 

6 

hi + 62 


4 


^^2 

4 


3^ ^2 

4 


\/ 2 c 2 


S- 


\Pic2 


s- 


3^2 „2 


<^2 <;,2 


(A3) 

(A4) 

(A5) 

(A6) 

(A7) 


and other field dependent masses given by. 


2 _ I 'J^hi ^2 /o 2 I a 2 \ I ^7,2 / A Q\ 

m<5 — —-I ^5^-1- —(S^ + A ) + ,(A8) 

" 1 ? = > (A9) 

, (AlO) 

= , (All) 


Field dependent masses matrix is given by. 


(Al) with G being goldstone fields. 


m\ s A = eigenvalues , 
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Appendix B: Scalar triple and quartic couplings 


gAAhi 


9AAh2 
gill hiii2 


gAAhihi 

gh2h2h2 


gAAh 2 h 2 

gh 2 h 2 hi 


gAAhih 2 


—2i{S2V COS (j )/4 + + d2xlA) simp) , 

—6i{Xv cos^ (p /4 + {d2x/4 + •\/2c2/12)) sin^ p 
+{82 x/4 + v^Ji/S) cos^ psinp 
+62x008^3111^^/4) , 

—2i{—S2Vsmp/4 + (•\/ 2 c 2/12 + d2x/4) cosp) , 
—2i{—3\v cos^ (/)sin p/4 + 3siiip^ cos0(d2a;/4 
+\p2c2/i-2) — { 62 X /2 + 'j26\/4) cos(/)sin^ p 
+ 62 v/ 2 siiipcos'^ p) , 

— 4 i{S 2 cos^ p/8 + d2 sin^ P/8) , 

—6i{—Xv sin^ p /4 + {d2x/4 + •\/2c2/12)) cos^ p 
+{ 62 x /4 + V^Ji/S) cosi^sin^ p 
—S 2 Vsmpcos'^p/ 4 ) , 

— 4 i {52 sin^ p/8 + d2 cos^ P/8) , 

— 2 i( 3 cos^ psinp{d2x/4 + •\/2c2/12) 

+3 sin^ p cos pXv /4 — {62X/2 + \f 26 \/ 4 ) 

X cos^(/)sin0— (52t'sin^(/)COS(/)/2) , 

—2i(s\VLpoosp{d2 —62)/^^ ■ (Bl) 


Here we briefly summarize the approach with notations 
adopted in Ref. |40] . With denoting the wall velocity, 
the plasma velocity u and the coefflcients ATi, Ki are de¬ 
fined by thermal averages. And the chemical potentials 
are expanded to second order in gradients. 

The source terms induced by the top quark are 

S^i = K-jd''TV ?, 

So = -VwKs{m^6')' +VwK^d'rn^{m^)' , 

Su = -Kiom^e'u[ . (Dl) 

in which the primes denotes derivatives with respect to 2 : 
(the direction perpendicular to the wall) and the CP-even 
first order perturbations ui, pi can be solved by 

Vi,Kip,[ -b ViiiK 2 {m^)'pi + u/ - r*°Vi = ViuKiiirri^y , 

-b VwK^u/ + VwKQ{wy)'ui + r*°‘ui = 0 . 

And, the Sg is the dominated one, effects induced by Su 
and are negligible. 

With the weak sphaleron rate |Sni, the strong 
sphaleron rate HQ], the top Yukawa rate [ 68 ], the top 
helicity flip rate, the Higgs number violating rate [55] . 
the W scattering rate, the quark diffusion constant m 
and the Higgs diffusion constant |72| being given by. 


Appendix C: Thermal masses used in EWPT 
analysis 


Thermal masses of bosonic field relevant to EWPT 
study are listed bellow, 

Hg = ( 3 gVl 6 + 5 VI6 + A/8 + y //4 + 62/24)T^ , 
^hh = ( 3 gVl 6 + <?'Vl6 + A/8 + y //4 + 62/24)T^ , 
Tlss = (^2/12 + S 2 / 48 )T^ , 

Haa = {d2/l2 + S2/48)T^ , 

H^ = n| = llg^QT^ 

H^ = Hi = 0 . (Cl) 


Appendix D: Transport equations 


P^, = 1.0 X 10-®T, 

Pj^ = 4.2 X 10"^r, 

^ m^{z,T) 
h , 


P^s = 4.9 X 10"^T, 

_ m‘l{z,T) 

m ggj, , 

■p -ptot 

^ h 5 

20 

Dh = ^ I (D2) 


the transport equations of chemical potentials, i.e., 
left-handed SU(2) doublet tops pt, 2 , left-handed SU(2) 
doublet bottoms fig, 2 , left-handed SU(2) singlet tops 
/.tt',, 2 (charge conjugation of tji), Higgs bosons Ph, 2 , are 
given by, 

2 + Suu, 2 + 3r^°*Uh^2 = 0 , 

~‘2Ki,hg'h,2 + 2Vii,K^ }iU/ 2 + 2r^°*Mft,,2 = 0 , 


For the case when the top quark acquires a com¬ 
plex mass during the electroweak phase transition, m = 
|m(z)|e*®(^\ the source terms and transport equations 
were derived in Ref. [ID] using the WKB approximation . 


2vii,Ki^}iy!h 2 + 2m/_2 ~ 8Ty(fj,t,2 + Mfe,2 + 2/rtc_2 + 2ph,2) 
—2Tiifih 2 = 0 , (D3) 
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iVwKi^ty^'t,2 + ^Vu]K2^t{lTl^yHt,2 + 3 M( 2 ~ ^^yifJ't,2 + /it‘=,2 + M/t,2) ~ 6rm(/Xt^2 + Mt‘=,2) ~ 3 rvt^(/Xt _2 ~ Mb,2) 
~3rss[(l + 9Ki^t)l-Lt,2 + (1 + 9Ki^i,)iii,,2 + (1 ~ 9-R'i_()/itc 2 ] = i , 

3Wu,-^l,bMb,2 + 3M'f, 2 — 3ry(/i;,_2 + Mt=,2 + Mb, 2 ) ~ 3rvi/(Mb,2 ~ Mi, 2 ) ~ 3rss[(l + 9Ki t)^t 2 
+ (1 + 9Ki^b)^b,2 + (1 ~ QiiTi t)/xtc 2 ] = 0 , 

3Wu)i^l,tMt':,2 + 3fu,i4r2,t(TOj )^/it<=,2 + 3?Xjc 2 ~ 3rj,(/i(_2 + Mb,2 + ‘^l^t‘^,2 + 2/X/i,2) ~ 6rm(/Xt_2 + Mt=,2) ) 

~3rss[(l + 9ifi t)/xt 2 + (1 + 9Ki b)Hb,2 + (1 ~ 9iiri_t)/itc_2] = tO'^rnlI , 

-3i4r4,tMi,2 + ivwKb,tuy2 + ivwk(,^t{'ml)'ut,2 + 3rt°*M4_2 = -3vn,Ks,t(mf9fy + ivwKg^t9'trnl{rnly 

-iKi^tkt-2 + Sv^kb^tuy2 + 3t'u;i^6,t(m4 )'utc^2 + 3r(°‘utc_2 = -iv^Ks^timfO'y}' + SvniKg^tO'tmtimt)' 

-3kw,tmXuyt^ ■ (D4) 
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